Suppression of dissipation in Nb thin films with triangular antidot arrays by random 

removal of pinning sites 
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The depinning current 7 C versus applied magnetic field B close to the transition temperature T c of 
Nb thin films with randomly diluted triangular arrays of antidots is investigated. Our experiments 
confirm essential features in I C (B) as predicted by Reichhardt and Olson Reichhardt [Phys. Rev. B 
76, 094512 (2007)]. We show that, by introducing disorder into periodic pinning arrays, I c can be 
enhanced. In particular, for arrays with fixed density n p of antidots, an increase in dilution Pd 
induces an increase in I c and decrease of the flux- flow voltage for B > B v = n p <&o. 

PACS numbers: 74.25.Qt, 74.25.Sv, 74.78.Na 
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The investigation of vortices in type-II superconduc- 
tors in the presence of tailored pinning potential land- 
scapes has attracted a lot of theoretical and experimental 
interest. On the one hand, vortices in superconductors 
may act as a model system in order to investigate general 
properties such as the dynamics and phase transitions in 
systems of interacting particles (e. g. colloidal suspen- 
sions P, [H, Wigner crystals [|| , charge density waves 
4 1 or various types of ratchets and Brownian motors 
E00])- O n the other hand, the ability to manipulate 
and control the static and dynamic properties of vortices 
is fundamental for superconducting device applications 

1- 

Modern lithography techniques allow the placement of 
artificial pinning sites into superconducting thin films 
with well-defined size, geometry and spatial arrange- 
ment. In case of periodic arrangements, enhanced vor- 
tex pinning was found for magnetic fields, at which the 
vortex lattice is commensurate with the pinning array 
@, EH, E3, [H, E3 ■ The enhanced pinning leads, e.g., 
to peaks in the critical depinning current I c at multiples 
of a so-called first matching field B p = n p $o; here the 
density of vortices carrying one flux quantum $0 = h/2e 
equals the density of pinning sites n p . However, at non- 
matching fields the vortex lattice is less pinned due to 
elastic deformations and formation of interstitial vortices. 
Hence, the question arises whether other arrangements - 
between the two extremes of periodic and random pin- 
ning arrangements - may lead to an enhanced vortex 
pinning over a broader range of applied magnetic field 
B. 

Recently, it has been shown by numerical simulations 
[lH UH and experimentally [13, E, Gil that a quasiperi- 
odic arrangement of pinnig sites produces additional 
commcnsurability effects and hence an enhanced pinning 
below the first matching field. A different proposal was 
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made very recently by Reichhardt and Olson Reichhardt 
pol ]. By molecular dynamics simulations they investi- 
gated periodic pinning arrays that have been diluted, by 
randomly removing the fraction Pd of pins, while keeping 
the pin density n p fixed. Such arrays are very interest- 
ing, since with increasing dilution the pinning potential 
undergoes a gradual transition from periodic to purely 
random. Therefore, this model is suitable to explore the 
intermediate region between order and disorder, as it is 
usually found in real world. Interestingly, the simulations 
showed that the introduction of some disorder leads to an 
enhanced critical current above the first matching field. 
In periodic pinning arrays the vortices sitting at the pin- 
ning sites form easy flow channels for interstitial vortices 
[2lj . while for randomly diluted pinning arrays channel- 
ing should be suppressed [2(j. This approach is also in- 
teresting from a general point of view, as the presence of 
disorder in competition with periodic potentials is also 
investigated in many other physical systems, e.g. two- 
dimensional conductors [22| , Ising ferromagnets [23[ , and 
Josephson Junction arrays I24l . [2a ]. For related recent 
work on vortex phases see [26[ . 

In this work, we present results on the experimental in- 
vestigation of vortex pinning and flow in superconducting 
Nb thin films containing randomly diluted triangular ar- 
rays of submicron holes (antidots) as pinning sites. We 
studied I C (B) at variable temperature T close to the su- 
perconducting transition temperature T c , and we com- 
pare pinning arrays with different dilution, considering 
two different scenarios: (i) "Scaled lattices": For differ- 
ent values of Pd, we fix the density n p of pinning sites. 
Accordingly, the lattice parameter (smallest separation 
between pinning sites) scales as a(Pd) = a(0)\/l — Pd- In 
this scenario a controlled transition from periodic to ran- 
dom arrangement of antidots is investigated, (ii) "Fixed 
lattices ": Here, we fix the lattice parameter a for differ- 
ent values of Pd- Accordingly, the density of pinning sites 
scales as n p (Pd) = (1 — Pd)n p (0). Here, with increasing 
Pd a transition to plain films (no antidots) is treated. 

Our experimental results confirm essential features as 
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predicted in [2(| : The fixed and scaled lattices show two 
different kind of matching effects, which differently de- 
pend on temperature. Furthermore, the scaled lattices 
show an enhancement of I c at magnetic fields above B p 
with increasing dilution Pd- This effect is caused by sup- 
pression of vortex channeling and can be also observed 
in the dynamic regime, i.e. by measuring current-voltage 
(IV) characteristics. 

The experiments were carried out on d = 60 nm thick 
Nb films which were deposited by dc magnetron sputter- 
ing in the same run on four separate Si substrates with 
1 /im thick SiC>2 on top. Patterning was performed by 
e-beam lithography and lift-off to produce Nb bridges of 
width W = 200 /mi and length L = 640 um. The bridges 
contain circular antidots (diameter D — 260 . . . 550 nm) , 
arranged in a triangular lattice that has been randomly 
diluted, with dilutions Pd=0 ("undiluted array"), 0.2, 
0.4, 0.6, 0.8, and 1 ("plain" film, without antidots). 
Each chip (#1 to #4) contains two or three sets (A, 
B, C) of bridges. Each set has six bridges with dif- 
ferent values for P c i = ... 1 . The antidot diameter 
D is kept constant within each set and varies from set 
to set. The chips #1 and #2 contain sets of bridges 
with fixed lattice parameter a = 1.5 /im. I. e. the den- 
sity of vertices of the corresponding triangular lattice is 
ni = -^^r = 0.5 fim~ 2 , which corresponds the "lattice 
matching field" B>i = ni$o = 1.1 mT (denoted as in 
[20j). For those two chips, the antidot density n p de- 
creases from 0.5 to 0.1 um~ 2 with increasing Pd from 
to 0.8. The chips #3 and #4 contain sets of bridges 
with scaled lattice parameters a(Pd) = 3.4 /im to 1.5 jum 
for Pd = to Pd = 0.8, respectively, in order to have 
a fixed antidot density n p = 0.1 /im -2 and "pin den- 
sity" matching field B p = 0.21 mT (denoted as B$ in 
[2g|). with N p ?s 12,500 antidots in each bridge. Below 
we present results obtained on bridges from sets #1-B 
(D = 300 nm), #3-B (D = 450 nm), #4-B (D = 360 nm) 
and #4- A (D = 260 nm). 

To characterize our devices, we first measured resis- 
tance R vs. T at B = and determined T c and normal re- 
sistance R n = R(T = 10 K) (with bias current 1 = 2 //A) 
of the different bridges on each chip. Due to the strong in- 
fluence of the reduced temperature t = T/T c (for T close 
to T c ) on the characteristic length scales, i.e. the Lon- 
don penetration depth X(t) and coherence length 
and on I c (t), the determination of T c plays an impor- 
tant role for the comparison and interpretation of the 
performance of pinning arrays with different Pd- We de- 
fined T c by linear extrapolation of the R(T) curves in 
the transition region to R = 0, i. e. T c marks the onset 
of resistance. For all samples we find T c w 8.5 K with 
a variation of a few mK within each set of bridges, and 
Rn = 5.0 to 5.817, depending on n p and D. Within 
the sets (D = const) with n p = const (scaled lattices), 
Rn varies by less than ±2 % from bridge to bridge. The 
plain film (from set #4-A) has R n = 5.0 £7, which yields 
a normal resistivity p n = R n dW/L = 9.4 /ificm. With 
the relation pi = 3.72 x 10~ 6 pflcm 2 [27| we estimate 
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FIG. 1: (color online) I C /I C (Q) vs B/B p for randomly diluted 
arrays at different reduced temperatures (I c increases with 
decreasing t). Comparison of the dependence on the antidot 
density n v for samples with Pd = 0.2 (a),(b) and Pd = 0.4 
(c), (d). 



for the mean free path £ = 4.0 nm. All I c values were 
determined with a voltage criterion V c = 1 /tV. 

Fig. Q] shows I C (B/B P ) patterns of four randomly di- 
luted bridges with similar antidot size at various temper- 
atures t = T/T c = 0.9995 . . . 0.9965. I c was normalized 
to its maximum value at B = 0. Two bridges (from 
set #4- A) have the same antidot density n p = 0.1 um~ 2 
and dilution P d = 0.2 (a) and P d = 0.4 (c). The two 
other bridges (from set #1-B) have larger antidot den- 
sity n p = 0.4 fjja with Pd = 0.2 (b) and n p = 0.3 /<m -2 
with P d = 0.4 (d). 

For the highest t = 0.9995, a clear peak in I C (B) in- 
dicates matching of the vortex lattice with the pinning 
array, as shown in Fig. [TJb). The position of the peak is 
located between B p and Bi. This indicates that indeed 
"pin density matching" is observed; however, as pointed 
out in [20] , for a diluted periodic pinning array at B = B p 
the vortex configuration contains numerous topological 
defects. Hence, commensurability effects at B p can only 
be observed when pinning is so strong, that the lattice 
distortion energy, associated with the deviation from an 
ideal triangular vortex lattice, can be overcome. This is 
most likely to be observed for the samples with higher 
density of pinning sites. Accordingly, for a given Pd, the 
matching peak is more pronounced in the samples with 
three and four times larger n p . 

With decreasing n p and t and with increasing Pd the 
peak in I C (B) gradually transforms into a shoulder- like 
structure, located close to B p . Following the evolu- 
tion of I C (B) with decreasing temperature shows that 
I C (B < B p )/I c (0) increases most whereas I C (B[) / I c (0) is 
almost independent of t. For the larger Pd = 0.4, this 
leads gradually to a transformation of the shouldcr-likc 
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FIG. 2: (color online) I c vs B/Bi at t = 0.9990 (a) and 0.9965 
(b) for arrays with fixed lattice parameter a and variable di- 
lution Pd. I C {B > Bi) decreases with increasing Pd- 



structure into a triangular-shaped I C (B) pattern with- 
out any indication of matching effects, either at B p or 
B~i. Nevertheless, I c for Pd — 0.4 is still significantly en- 
hanced over I c for samples without antidots, as will be 
shown below. 

In the following, we directly compare pinning arrays 
with different dilution Pd at the same reduced temper- 
ature. Fig. [2] shows I C (B) patterns of samples with 
fixed lattice parameter a — 1.5 /xm (n; = 0.5 /im -2 ) at 
t = 0.9990 (a) and t = 0.9965 (b). The B-axis is normal- 
ized to the lattice matching field B{ , which is the same for 
all perforated bridges within this set. The sample with 
Pd = shows pronounced peaks in I C {B) which are lo- 
cated at ±B; and ±25;, indicating a saturation number 
n s > 2 [28l . [29| for both temperatures. In the samples 
with small dilution (Pd = 0.2 and 0.4) we also find peaks 
in I C {B) for the higher temperature shown in FigJ^a). 
The matching peaks are significantly broader than the 
matching peaks of the undiluted bridge, and they are lo- 
cated at magnetic fields between Bi and B p (Pd). This is 
also visible in Fig. [TJb) and (d). With increasing Pd 
we find a gradual transition of the I C (B) patterns at 
B < Bi from the undiluted array (Pd = 0) to the plain 
film (P d = 1). Interestingly, for t = 0.9990 [c.f. Fig^a)] 
the diluted sample with Pd = 0.2 shows a higher I c than 
the undiluted sample (Pd = 0) for B < Bi, and very sim- 
ilar I c values for all other fields, except for the matching 
fields Bi and 2B t . 

Fig. [3] shows I C (B/B p ) patterns of bridges with dif- 
ferent P d = . . . 1 at t = 0.9990 (a) and * = 0.9945 
(b). In Fig. Ola) the data from set #4- A (D = 260 nm) 
are complemented by results from two undiluted bridges 
with D = 360 nm (#4-B) and D = 450 nm (#3-B), in or- 
der to demonstrate the effect of antidot size. All bridges 
with Pd — and different D show qualitatively the same 
I C (B) patterns. The major difference is a slight increase 
in I c at B < B p with increasing D. This dependence can 
be explained with the increase of the pinning strength 
of the antidots with increasing D. For B < B p , each 
vortex can be captured by an antidot, and hence the in- 
creasing pinning strength with D leads to an increase of 
I c . In contrast, for B > B p , vortices occupy interstitial 



FIG. 3: (color online) I c vs B/B p at t = 0.9990 (a) and 0.9945 
(b) for arrays with scaled lattice parameters a(Pd) (fixed anti- 
dot density n p = 0.1 /im~ 2 and matching field B p ): Compar- 
ison of samples with different dilutions Pd and with different 
antidot diameters D for Pd = 0. I C (B~B P ) decreases with 
decreasing D (for Pd = 0) and with increasing Pd- Vertical 
solid lines indicate ±Bi(Pd)/B p , labeled with corresponding 
values for Pd- 



pinning sites where they are weakly pinned. Hence, I c 
is drastically reduced and determined by the motion of 
interstitial vortices, which should not depend on the an- 
tidot size. This is exactly what we find experimentally, 
i.e. I C (B > B p ) is independent of D for the samples with 
P d = 0- 

With increasing Pd (decreasing a in order to keep the 
antidot density n p constant) the shape of the I C (B/B p ) 
patterns is strongly affected. The diluted arrays show 
"lattice matching" effects around It is interesting to 
note that with decreasing t the matching effects become 
less pronounced, and the I C (B) pattern approaches more 
and more a triangular shape. This shape is reminiscent 
of the I C (B) pattern of randomly arranged antidots [l8|. 
This can be understood, as with increasing Pd and fixed 
n p , the antidot arrangement approaches that of a random 
arrangement. 

For Pd = from set #4-A we have no data for the full 
range of t. as the bridge was damaged after taking first 
data. Hence, to facilitate the comparison of different P d 
at t = 0.9945 [Fig.^b)], we also show data from another 
bridge with Pd = (from set #4-B). For both tempera- 
tures (t = 0.9990 and t = 0.9945) we find a decrease in 
I c with increasing Pd for B < B p . However, for B > B p 
the diluted arrays show an enhanced I c as compared to 
the undiluted samplc(s). Both observations are in qual- 
itative agreement with the simulations in Ref. [20| . We 
do find that the enhancement of I c above B p persists up 
to the highest dilution Pd = 0.8. This enhancement was 
explained in po| with the suppression of channeling of 
interstitial vortices for fields in the range Bi < B < B p . 
For undiluted arrays, this channeling effect causes the 
rapid decrease of I c with increasing B slightly above B p 
[c. f. Fig. [5]. Our experimental data clearly confirm that 
this rapid drop in I c at B p is absent for the diluted arrays. 

The suppression of channeling should also be visible 
in the current-voltage-characteristics. Fig. 2] shows the 
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FIG. 4: (color online) V(I) curves at B = 2B P and t = 0.9945 
for arrays with scaled lattice parameters a(Pd) (fixed antidot 
density n p — 0.1 fim~ 2 and matching field B p ). Upper inset 
shows magnification of V vs I for Pd = with the critical 
depinning current I c ,\(Ic,2) of interstitial (pinned) vortices. 
Lower inset shows magnification of V vs I at small voltages; 
V(I >1 mA) decreases with increasing Pd, except for the plain 
film (Pd = 1) with lowest 7 C . 

V(/)-curves of differently diluted samples at B = 2B p 
and t — 0.9945. The V r (/)-curves directly correspond 
to the I C (B) pattern shown in Fig. [3fb). For the undi- 
luted bridge at B = 2B p one expects the same number 
of vortices sitting in the antidots and at interstitial posi- 
tions. In this case we find two critical depinning currents 
Ic,ijc,2 in the V^(/)-curves [c.f. upper left inset of Fig. [4]. 
I Ci i corresponds to the current above which a finite volt- 
age appears. This voltage is caused by the motion of 
weakly pinned intcrstitials. At a higher current / = Z Cj 2 
the slope of the TV-curve changes. This is due to the 
depinning of vortices sitting in the antidots. The value 
of I c ,2 can also be found in the I C (B) pattern shown in 
in Fig. [3] (b). fits quite well to the critical current 
at the first matching field I C (B P ). All diluted samples 
(Pd = 0.2 . . . 0.8) have a similar I c i but do not show a 
pronounced slope change at I c .2- The lower right inset in 
Fig. [4] clearly shows that with increasing Pd the voltage 



due to flux motion decreases. I. c.. in the diluted pinning 
arrays dissipation is reduced, due to the more effective 
suppression of vortex channeling. 

In conclusion, we experimentally investigated Nb thin 
films with triangular arrays of antidots. which have been 
randomly diluted, by measurements of the critical cur- 
rent I c vs. applied magnetic field B and current-voltage 
(IV) characteristics close to the transition temperature 
T c . The antidot lattices could be tuned to find two differ- 
ent matching effects, related to the antidot density and 
to the lattice parameter of the antidot lattice, as pre- 
dicted in [2(j. For samples with fixed lattice constant, 
with increasing dilution Pd a gradual transition from a 
periodic pinning array to a plain film without pinning 
sites has been observed. Obviously, with increasing Pd 
the critical current decreases. However, very close to T c , 
for small dilutions (Pd = 0.2) we do find a broad peak in 
I C (B) located between B p and Bi, corresponding to an 
increase in I c by removing 20 % of the pinning sites. We 
speculate that this counterintuitive effect is due to the in- 
troduction of disorder; an understanding of this effect is 
still lacking, and deserves further investigations. On the 
other hand, for samples with fixed antidot density, an 
increasing dilution corresponds to a gradual transition 
from a periodic to purely random distribution of pinning 
sites. Our experiments clearly show an enhancement of 
I c for magnetic fields above B p with increasing Pd- This 
was the main prediction in Ref. [2fj| and can be explained 
with the suppression of channeling of interstitial vortices. 
This effect is also observed in TV-measurements, i.e., the 
suppression of channeling causes an increasing reduction 
in the flux- flow voltage with increasing Pd- As a conse- 
quence, the concept of introducing disorder by randomly 
removing pinning sites in tailored periodic pinning arrays 
seems to provide a feasible way for enhancing the critical 
current in superconductors for magnetic fields above the 
matching field B p . 
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